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Abstract 

For a class of one-dimensional Schrodinger operators with polyno- 
mial potentials that includes Hermitian and PT-symmetric operators, 
we show that the zeros of scaled eigenfunctions have a limit distribu- 
tion in the complex plane as the eigenvalues tend to infinity. This limit 
distribution depends only on the degree of the polynomial potential 
and on the boundary conditions. 
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1. Introduction 

We begin with an eigenvalue problem of the form 

- y" + P(x)y = Ay, y(-oo) = y(oo) = 0, (1) 

where P(x) = x d + . . . is a real monic polynomial of even degree d. The 
boundary condition is equivalent to y G L 2 (R). It is well-known that the 
spectrum of this problem is discrete, all eigenvalues are real and simple, and 
they can be arranged into an increasing sequence Ao < Ai < . . . — > +oo. 
Moreover, 
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where B and V are the Euler's functions. A general reference for these facts 

isnu. 

Eigenf unctions y n are entire functions of order (<i+2)/2, and in this paper 
we study the distribution of their zeros in the complex plane when n is large. 
When d = 2 and P is even, we have y n (x) = H n (x) exp(— x 2 /2), where H n 
are the Hermite polynomials, and asymptotic distribution of zeros is known 
in this case in great detail [IB] . The case d = 4, which is called sometimes 
an anharmonic oscillator or a double well potential, was studied much, but 
most attention was payed to the properties of eigenvalues, rather than the 
properties of eigenfunctions. 

In [3J, we proved that for d — 4 and even P, all zeros of all eigenfunctions 
belong to the union of the real and imaginary axis. The main result of the 
present paper implies that after an appropriate rescaling of the independent 
variable, zeros of eigenfunctions have a limit distribution in the plane, which 
depends only on d. This will be derived from the asymptotics of log \ Y n \ as 
n —>■ oo, where Y n is a properly rescaled eigenfunction. 

Let us consider real zeros first. According to the Sturm-Liouville theory, 
y n has exactly n real zeros, all of them simple. A classical argument shows 
that all these real zeros belong to the interval (a~,a+) where a~ and a+ 
are the smallest and the largest roots of the equation P(x) — A n = 0. The 
asymptotics 

le^l ~ X l/d 

suggests to define the rescaling: 

Y n (z) = y{\]! d z). (2) 

Now we denote by v n the counting measure of the roots of Y n (for every 
set X C C, v n {X) is the number of roots on Y n in E). Our main result, 
Theorem 2 below, has the following corollary: 

— — > Cd\l 1 — x d dx, — 1 < x < 1, 

n 

where the convergence is the usual weak convergence of measures: /i n — > 
fi means that / <pd^ n — > J <pdfi for each continuous function with compact 
support. The normalizing constant is 

c rf = ~fl(3/2,l/d). 
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For d = 2 the limit distribution is the "semi-circle law", the well-known 
asymptotic distribution of zeros of Hermite's polynomials. 

A theorem of Hille [9J Theorem 11.3.3] implies that for every r 6 (0, 1) 
there exists no(V) > such that for all n > no, all zeros of Y n in the disc 
{z : \z\ < r} are real. 

Zelditch [IB] studies Laplace-Beltrami eigenfunctions on real analytic 
manifolds. Under certain conditions on the manifold, he extends the eigen- 
functions to a complex neighborhood of the manifold and obtains a limit 
distribution of their zeros. Asymptotics in the complex domain help to study 
the distribution of real zeros. The same happens in our case. 

Passing to the consideration of complex zeros, we make the same rescaling 
([2]), and consider the rescaled counting measures 

/i« = v n jn (3) 

of zeros of Y n in the complex plane. 

Our main result will imply that these measures /i converge weakly to an 
explicitly described limit measure, which depends only on d. 

Our results also apply to the so-called PT-symmetric Schrodinger op- 
erators which were intensively studied in the recent years [21 [13]. Let P 
be a complex polynomial of degree d with the property P(—z) = P(z). 
Schrodinger operators with such potential P are called PT-symmetric. Ev- 
ery PT-symmetric potential can be written as P{z) = Pi(iz), where Pi is a 
polynomial with real coefficients. A real potential P is PT-symmetric if and 
only if P is even. 

Following K. Shin [13J, we consider the generalized eigenvalue problem 
which contains both self-adjoint and PT-symmetric cases. 

-y" + Py = \y, (4) 

with 

P{z) = {-l)\iz) d + Y,a k z\ 

k=l 

where the coefficients are arbitrary complex numbers, and the boundary 
condition is 

y(r C ")-0, r^oo for 9 = -| ± (5) 
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where 1 < £ < d — 1. The self-adjoint problem (pQ) corresponds to the case 
that d is even, i = d/2, and are real. The usual boundary condition 
imposed on a PT-symmetric potential corresponds to £ — 1 PQI2] ■ 
The following result belongs to Sibuya [14J and K. Shin [13J. 



Theorem A. The spectrum of the problem §4§, Q is discrete, and 
to each eigenvalue corresponds a one- dimensional space of eigen] 'unctions. 
Eigenvalues satisfy 

( y/HX3/2 + l/d)n (R . 

An ~ U(W<*)r(i + i/d) J ' n ^°°" (6) 

So we see from the asymptotics that the eigenvalues are "asymptotically real" 
in the sense that their arguments tend to zero. Shin proved that in the PT- 
symmetric case (that is when at = {ibk) k with real bk), all eigenvalues but 
finitely many are actually real. 

Our main result, will imply that zeros of eigenfunctions y n of the problem 
OH), ©, when properly rescaled as in (jSJ) have a limit distribution that de- 
pends only on d and I. The support of the limit distribution consists of some 
Stokes lines (which we later define precisely), and our result can be consid- 
ered as a rigorous confirmation of the results of numerical computations of 
Bender, Boettcher and Savage [TJ. 

Functions Y n satisfy the differential equations of the form 

Yn = K /d (PM d z) - K)Y n = k 2 n ((-l) e (iz) d - 1 + o(l))Y n , (7) 

where k n = A^ +2 ^ ( - 2d - ) — > oo, as n — > oo. Here we choose the branch of 
\( d + 2 )/( 2d ) which is positive on the positive ray. 

It is well-known (see, for example [HE]) that the asymptotic behavior, as 
n — > oo, of solutions of such differential equations depends on the quadratic 
differential 

Q d/ (z)dz 2 = ((-lY(tz) d - l)dz 2 . (8) 

We recall the terminology and some known facts about such quadratic differ- 
entials. For the general theory of quadratic differentials we refer to [TH [T5] , 
and for applications to differential equations to [HE]. 

Let Q be an arbitrary polynomial. The zeros of Q are called the turning 
points. The curves on which Q(z)dz 2 < are called the (vertical) trajectories. 
Each branch of the integral / \fQdz maps trajectories into vertical lines. The 
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trajectories make a foliation of the plane with the singularities at the turning 
points. The leaves of this foliation are maximal smooth open curves on which 
Q(z)dz 2 < holds. Each leaf is homeomorphic to an open interval, and each 
end of a leaf is either at infinity or at a turning point. Those leaves which 
have at least one end at a turning point are called the Stokes lines. A Stokes 
line whose both ends are turning points is called short. At each simple zero 
of Q exactly three Stokes lines converge, and they make equal angles of 2tt/3 
between them. 

The components of the complement of the union of turning points and 
Stokes lines will be called the Stokes regions. These regions are unbounded 
and simply connected. The closure of each Stokes region (in C) is homeo- 
morphic either to a closed half-plane or to a closed strip. We say that these 
regions are of the half-plane type or of the strip type, respectively. 

Turning points, Stokes lines and Stokes regions make a decomposition of 
the plane which we call the Stokes complex. 

Our first result is the topological description of the Stokes complex cor- 
responding to the differential (jSj). 

Theorem 1. The Stokes complex of Qd/dz 2 is symmetric with respect to the 
imaginary axis and has the following property: Every turning point v which 
does not lie on the imaginary axis is connected by a short Stokes line with 
the turning point —v symmetric to v with respect to the imaginary axis. 

It is easy to see that this theorem yields a complete topological description 
of the Stokes complex. The turning points are the roots of the equation 
(— l) e (iz) d = 1. All these roots are simple, so three Stokes lines meet at each 
turning point. 

If v is a turning point in the (open) right half-plane, then there is one 
short Stokes line from v to —v. The other two Stokes lines originating at v go 
to infinity. Indeed, if one of those two were short, we would obtain a bounded 
Stokes region which is impossible. These two Stokes lines are contained in the 
right half-plane (otherwise they would intersect the symmetric Stokes lines 
originating at — v which is impossible.) The picture in the left half-plane is 
symmetric to the picture in the right half-plane. 

If v belongs to the imaginary axis, that is v = ±z, then all three Stokes 
lines originating from v are unbounded. One of them coincides with a ray of 
the imaginary axis (±i, ±ioo) and the other two are symmetric with respect 
to the imaginary axis. 

See Figures 1-9 where the Stokes complexes of Qd/dz 2 for d — 2, 3, 4 and 
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6 and various values of i are represented. 
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Fig. 3. d = 4, I = 1. 





Fig. 6. d = 6, I = 1. 




Fig. 7. d = 6, £ = 2. 




Fig. 8. d = 6, £ = 3. 
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Fig. 9. d = 6, t = 4. 



It is easy to see that each Stokes complex has exactly d + 2 Stokes regions 
of the half-plane type. These regions are asymptotic to the sectors which are 
bounded by the Stokes directions 

{e ie : Re jH^ y/(-iy(uydt = 0} that is g = -| + g£L±|Q . 

The bisectors of the angles between adjacent Stokes directions are called the 
anti-Stokes directions. Each of the d+2 anti-Stokes directions approximately 
bisects some Stokes region of the half-plane type. Thus the boundary con- 
dition (jSJ) indicates that y(z) — > as z — > oo on two anti-Stokes directions 
symmetric with respect to the imaginary axis. This boundary condition sin- 
gles out two Stokes regions of the type of half-plane: the right region u + and 
the left region uj~, so that the eigenfunction tends to zero along the bisectors 
of oj~ and u + . In Figures 1-9 the regions omega + and u~ are marked by 
+ and — signs, respectively. Figures 1, 4 and 8 correspond to self-adjoint 
problems. 

It follows from the topological description of the Stokes complex of Qd,tdz 2 
above that the closures of u + and u~ are disjoint. 

To state our main result, we need to define certain Stokes lines in the 
Stokes complex of Qd/dz 2 which we call exceptional. There will be exactly 
one exceptional Stokes line originating at each turning point. 

Let v + and v~ be the turning points on the boundaries of uj + and u~ , 
respectively. Then v + and v~ do not belong to the imaginary axis, so The- 
orem 1 implies that there is a short Stokes line Eq = (v~,v + ). This short 
Stokes line will be exceptional. For example, for the self-adjoint problem ([T]) 
we have E = (—1, 1). 
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If v = % or v = — % is a turning point, then the Stokes line (i, +ooi) or 
(— i, —ooi) is exceptional. 

For the rest of turning points, exceptional Stokes lines are defined as 
follows. The complement of the set u + U oj~ U Eq, where the bar stands 
for the closure, consists of two components, we call these components D + 
(containing the positive imaginary ray) and D~ (containing the negative 
imaginary ray). Let v be a turning point in D + and not on the imaginary 
axis. Let V and L" be the two unbounded Stokes lines originating at v. Of 
these two Stokes lines we choose that one which lies between the other one 
and the positive imaginary axis, and call this chosen line exceptional. (Any 
family of disjoint curves tending to infinity in D + can be linearly ordered, 
for example anticlockwise; we used the word "between" in the sense of this 
order). Similarly, if v is a turning point in D~ and not on the imaginary axis, 
then of the two unbounded Stokes lines originating from v, the exceptional 
one is that which lies between the other one and the negative imaginary axis. 

Exceptional Stokes lines are shown as bold lines in our figures. 

Let E be the union of all exceptional Stokes lines and all turning points. 
We call E the exceptional set and denote 

Q = C\E. 



Then f2 is a doubly connected region, and y (— iy(iz) d — 1 has two single- 
valued holomorphic branches in Q. Consider the harmonic function in Q 

u{z) = Re J* ^(-iy(U) d -ldt = Re J* ^Ch/t)dt. (9) 



Here we choose the branch of the yQd,e i n such a way that u(z) — > — oo as 
z — > oo along the anti-Stokes directions in uj + and u~ . (The normalization 
(J5j) and the top coefficient of the potential P in (Tj0) are chosen to make the 



choice of such branch of JQd,£ possible.) 

The integral in the definition of u has one period corresponding to loops 
in Q around the short Stokes line Eq = (v~,v + ), but this period is pure 
imaginary, because Qd/dz 2 < on E Q , thus the real part of this integral is a 
well defined harmonic function in Q. In fact u is continuous and subharmonic 
in the whole plane (that the jumps of the integral on the exceptional Stokes 
lines is pure imaginary can be seen from the very definition of the Stokes 
lines). Now we can state our main result. 
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Theorem 2. Lety n be an eigenfunction of the problem Q, Q corresponding 
to the eigenvalue \ n , and normalized so that \y n (0)\ + \y' n (0)\ = 1. PutY n [z) = 
y n (\ l J d z). Then 

- log \Y n \ — > c d m(z), n — > oo, 

uniformly on compact subsets of Q, and also in the sense of Schwartz distri- 
butions D' in C. Here u is the function defined in and 

v/?r(3/2 + l/d) 



sin(&r/d)r(l + l/d)" 



Uniform convergence implies that the functions Y n have no zeros on a 
compact C fi when n > tiq{K). This can be made more precise by 
replacing a compact set by an unbounded subset of Q which we describe 
below as an "admissible set" . Each compact subset of Q is contained is some 
admissible set. 

The Laplace operator is continuous in D', and the Riesz measure fi = 
(27t) _1 Am can be easily explicitly computed: 

Corollary. The normalized counting measures \x n = v n jn of the zeros ofY n 
converge weakly to the measure 

A* = c d ,e\J\QdA z )\dz\, 

supported on E. 

In particular, in the self-adjoint case ([1]) the limit density on (—1, 1) 
is given by the formula qVT advertised. Paper pQ contains nice 

pictures showing the zeros of the rescaled eigenfunctions of PT-symmetric 
operators clustering around E . 

2. Proofs 

Proof of Theorem 1. It is convenient to make the change of the indepen- 
dent variable z > iz which reduces our quadratic differential to the form 
(±z d + l)dz 2 . It is clear that the Stokes complex of this new differential is 
symmetric with respect to the real line. All we have to prove is that there is 
a short Stokes line from every turning point in the upper half-plane to the 
complex conjugate turning point. The proof is performed in several steps. 
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Step 1. Let P(z) = (z d + 1). We prove that the Stokes complex of 
P[z)dz 2 contains a short Stokes line between the turning points v± = e ln ^ d 
and V-\ = vl = e~ ln ^ d . To find the directions of the Stokes lines originating 
at vi, we write P(z) = dvf~ 1 (z — V\) + o(z — i>i), and 

r jp{t)dt = -Vdv{ d - i)/2 ( Z - Vl f 2 + o{ Z - Vl f 2 . 

Jvi 3 

So on the tangents to the Stokes lines meeting at V\ we have (z— v^vf^ 1 < 0. 
Since vf^ 1 = — l/t>i, this is equivalent to 

(z — f i) 3 G t>iR + . (10) 

Consider the sector S bounded by the segments [0, 1], [0,t>i] and the arc of 
the unit circle [1, Vi]. Our polynomial P maps S into the first quadrant, thus 
for z G S, we have arg(^ d + 1) G (0, 7r/2). This means that the direction 
of the vertical foliation P(z)dz 2 < in S satisfies 

tt/4 < (p + nk < tt/2. (11) 

There is one Stokes line originating at v± with the direction in this interval. 
According to (fit)]) , its direction at v± is 7r/(3<2) — 27r/3. Because of the lim- 
itations (ITTj) on the direction of this line in 5", it can leave 5 only through 
the interval (0, 1). On this interval it has to meet the symmetric Stokes line 
from v^i. Thus V\ and t>_i are connected with a short Stokes line. 

Step 2. The vertical foliation of the differential Pg(z)dz 2 = e~ 2e P(e~ e z)dz 2 
is obtained from the vertical foliation of P(z)dz 2 by counterclockwise rota- 
tion by 9. In particular, the Stokes complex of Pg(z)dz 2 contains a short 
Stokes line between the adjacent turning points e l ^ e±n ^ d \ 

Step 3. We prove that every turning point Vk = e m ( 2k ~ 1 ^ d of P{z)dz 2 in 
the (open) first quadrant is connected to the symmetric turning point by a 
short Stokes line. This we prove by induction in k. For k = 1 the statement 
was proved in Step 1. Suppose it holds for k < m — 1, and let L m _i be the 
Stokes line from f m _i to its conjugate. 

The tangents to the Stokes lines originating at v m satisfy the equation 
(z — v m ) 3 G f m R+ similar to ffTOl) . In particular there is one Stokes line 
L m starting at v m in the direction (p m = 7r(2m + l)/(3d) — 27r/3, so that 

-27T/3 <(p m < -TT/2. 
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From Step 2 applied with 9 = 2nm/d, there is a Stokes line U m for the 
differential P e dz 2 connecting v m and i> m -i- Its tangent at v m has direction 

i/j m = 0! + 9 = n/d - 2tt/3 + 2nm/d = n(2m + l)/d- 2tt/3 > <p m . 

We are going to show that L m never intersects U m in the open unit disc. 
Proving this by contradiction, we suppose that z is a point of intersection, 
and denote by L' rn and U' m the pieces of L m and U m from v m to z. Then the 
integrals 

f ^P{z)dz and f yjp e (z)dz 

are both non-zero (because the imaginary part of such integral is monotone 
along a Stokes line) and both pure imaginary. But P e = e - 2nirn / d P for our 
choice of 9, so we conclude that 

/ yfpjzjdz ± f Jpjzjdz, 

as these integrals have different arguments and are non-zero. This is impos- 
sible because the integral of y/Pdz should be zero over any closed curve in 
the unit disc. Thus L m and U rn do not intersect in the open unit disc. 

For all z in the unit disc, the direction of the vertical foliation of P(z)dz 2 
belongs to the interval n/4 < <fi + irk < 37r/4. Hence L m cannot leave the 
upper half of the unit disc through the arc [v m , —1] of the unit circle. 

Moreover, L m cannot intersect the Stokes line L m _i which exists by the 
induction assumption. 

This leaves only one possibility that L m leaves the upper half of the unit 
disc through the real axis. Then, by symmetry, it should continue to f_ m . 

Step 4- Now we prove the same for the differential P_(z)dz 2 = (—z d + 
l)dz 2 , namely that every turning point in the first quadrant is connected by 
a short Stokes line with its conjugate. First we prove that the Stokes line 
Li starting at the turning point Z\ = e 2m l d in the direction 27r/(3rf) — 27r/3 
cannot intersect the Stokes line of P^/ddz 2 = e~ 2n / d P^(z)dz 2 connecting z\ 
with 0, hence the only possibility for L\ to leave the upper half of the unit 
disc is through the real line, so it should proceed to z~[ by the symmetry. 
Then we proceed by induction on m to prove that there are Stokes lines 
for P_dz 2 connecting the turning points z m in the first quadrant with their 
complex conjugates. The induction step is similar to the one in Step 3 and 
we leave it to the reader. 
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Step 5. That every turning point in the left half-plane is connected by a 
short Stokes line with its complex conjugate is proved by the change of the 
independent variable z i— > — z. 

This completes the proof of Theorem 1. 

Proceeding to the proof of Theorem 2 we begin with some rough esti- 
mates. 

Lemma 1. Consider a normalized solution of the differential equation 

y" = h 2 Qy, 2,(0) = y , y'(0) = y u 

where h > and Q is a holomorphic function satisfying \Q(z)\ < M for 
\z\ < R. Then we have 

\y(z)\ < max{|y |, |yi|}exp(/iAfi2). 

Proof. We put ui\ = hy and u>2 = y', to obtain a matrix equation 

w'(z) = A(z)w(z), 

where 

W= {Z) ^ A= {hQ o)- 

This implies 

iiw(*)n<iiw(o)n+ r \\A(t)\w\w(t)\\dt, 

Jo 

where we use the sup-norms 

|| w || = max{|wi|, \w 2 \}, 

and 

\\A\\ = max{|an| + |ai 2 |, |a 2 i| + |a 2 2 1 } < hM. 

Applying Gronwall's lemma P, Thm. 1.6.6], we obtain the result. 

Applying this result to the equation (J7j), and taking into account the 
asymptotics of the eigenvalues from Theorem A, we conclude that if the 
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sequences (1/n) log ^(zo)! and (1/n) log |Y^(z )| are bounded at some point 
z G C, then in every disc \z\ < R the family of subharmonic functions 

u n = -\og\Y n \ (12) 
n 

is uniformly bounded from above. By a well-known theorem (see, for ex- 
ample, [TO, Theorem 4.1.9]) it follows that from every sequence of these 
functions one can select a subsequence which either converges uniformly on 
every disc to — oo, or converges in D'(C) (or in L 1 1 oc (C) which is equivalent 
for such families of subharmonic functions). Convergence to — oo is excluded 
by normalization of our eigenf unctions. 

Thus the limit functions of our family are subharmonic, and the same the- 
orem PHI Theorem 4.1.9] says that the Riesz measures of u n weakly converge 
to the Riesz measures of u. 

Thus all we need to do is to identify the possible limit functions u, and 
to show that there is only one. This will be done by proving the uniform 
convergence part of Theorem 2. Notice that subharmonic functions are upper 
semi-continuous, so the restriction of u on Q defines u uniquely in the whole 
plane. 

To prove the uniform convergence in Theorem 2, we need a finite open 
covering of Q with the so-called canonical regions [I]. 

Consider a quadratic differential Q(z)dz 2 with arbitrary polynomial Q. 
The multi-valued function 

C(z) = f \[Q(t)dt (13) 

has holomorphic branches in each Stokes region. These branches are defined 
up to post-composition with a transformation z \— > ±z + c, where c is an 
arbitrary constant. Each such branch maps its Stokes region univalently 
onto some right or left half-plane, or onto a vertical strip, depending on the 
type of the region [U [5] . 

A region D, which is a union of Stokes regions and Stokes lines is called 
a canonical region if there is a holomorphic branch of ( in D that maps D 
onto the plane with finitely many slits along some vertical rays. 

It is clear that a canonical region contains exactly two distinct Stokes 
regions of the half-plane type. It is easy to state a criterion for a pair of 
Stokes regions of the half-plane type to belong to some canonical region: 
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Stokes regions Di and D 2 belong to a canonical region if and only if there 
is a curve from D\ to D 2 that does not pass through the turning points, 
intersects the Stokes lines transver sally, and intersects each component of 
the 1-skeleton of the Stokes complex at most once. 

For example, in Figure 6 there is a canonical region containing the Stokes 
regions —,5,6, another one containing —,5,4,1 and so on, but there is no 
canonical region containing — and +. 

It is easy to conclude from Theorem 1 that in the case that Q = Qd,e, for 
each Stokes region D, other than uj~, there is a canonical region containing 
u + and D. Similar statement holds with oj~ and u> + interchanged. Moreover, 
for Q = Qd,£, canonical regions containing u + cover Q, and canonical regions 
containing ui~ cover fl as well. So the statement about uniform convergence 
in Theorem 2 is enough to prove for all canonical regions containing u + or oj~ . 

Let Q(z)dz 2 be a quadratic differential with a polynomial Q, and let a 
number s > be given. A smooth curve 7 : (—1,1) — > C will be called 
s-admissible if it has the following properties: 

(i) 7(t) — * 00, as \t\ — > 1, 

(ii) For every t G (—1, 1), the (smaller) angle between ^'(t) and the vertical 
trajectory at ^(t) is at least s. 

(hi) For every t G (—1,1), the distance from j(t) to the set of turning points 
is at least s. 

A subset K C Q, will be called s-admissible, if K is a union of s-admissible 
curves. 

A subset K C VL or a curve 7 in fl will be called simply admissible if it 
is s-admissible for some s > 0. It is easy to see that interiors of admissible 
subsets of canonical regions containing oj + or oj~ cover Q. Two Stokes regions 
of the half-plane type belong to some canonical region if and only if there is 
an admissible curve passing through these two Stokes regions. 

Now we have to study the behavior of the Stokes complex under small 
perturbations of the quadratic differential. For our purposes the following 
lemma will be sufficient. 

Lemma 2. Let Q(z, h) be a polynomial of degree d with respect to z whose 
coefficients are continuous functions of h, for h in a neighborhood of 0, and 
suppose that deg^ Q(z, 0) = d, and that K is an admissible subset of some 
canonical region of Q(z,0)dz 2 . Then there exists ho > 0, such that for 
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\h\ < h , the differential Q(z,h)dz 2 has a canonical region in which K is 
admissible. 



Proof. Let 7 be an admissible curve for Q(z, 0)dz 2 . This means that 7 
does not pass through the zeros of Q(z,0) and that the angles between the 
tangent lines to 7 and the vertical trajectories of Q(z,0)dz 2 are bounded 
from below by some s > 0. It is clear that for every compact piece of 7 
these conditions persist for the vertical trajectories of Q(z, h)dz 2 for small h. 
Thus one only has to consider a neighborhood of 00. Suppose that Q(z, h) = 
a(h)z d + . . .. Then there exist ho > and R > such that 

I argQ(z, h) — d&rgz — arga(0)| < s/2, 

whenever \z\ > R and \h\ < ho- The angle between 7' and the vertical 
direction is 

/ 7T argQ 
arg7 - 2 + — 2~ + 7r/c - 

so if this angle is at least s for h = 0, then it is at least s/2 for \h\ < h and 
I z I > R. This proves the lemma. 

The following lemma is essentially well-known, though we could not find 
a convenient reference in the literature, so we include a proof. This lemma 
constitutes the essence of the Liouville method, which is known to physicists 
as the WKB method (see [TJ [12] for the general discussion of the method, 
including history). The closest statement to ours is the one in [1], but our 
proof is different, and it is based on Liouville's transformation^ 

Lemma 3. Let Q be an arbitrary polynomial, and consider the differential 
equation 

-y" + h 2 Q(z)y = 0, (14) 

where h > is a parameter. Let D be a canonical region for Q(z)dz 2 , con- 
taining a Stokes region D\ of the half-plane type. Let ( = §(z) be a branch 
of 



'Q{w)dw 

which maps D onto the plane with vertical slits, and such that D\ corresponds 
to a left half-plane, and let K C D be an s-admissible set. Then, for each 



1 We do not share the opinion of Fedoryuk Ch. II, §1, sect. 6] that the Liouville's 
transformation is "less convenient" for complex functions than for the real ones. 
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h > h (s,Q), there exists a solution y^(z) of the differential equation fli^p , 
which satisfies 

y* h (z) = Q~ l /\z) exp(M>(^))(l + e(z, h)), (15) 

where 

\e(z,h)\<-\, zeK. (16) 
n — h 

Moreover, one can take 

5 Q' 2 Q" 



h {s,Q) = sup / 

J -V 



16 Q 3 4Q 2 



\d(\, (17) 



where the sup is over all s -admissible curves in K , and derivatives in l \17\ ) 
are with respect to z. 



Remark. It is easy to see that the integral in ffTTl) is absolutely convergent. 
Assuming that the leading coefficient of Q has absolute value 1 and all roots of 
Q belong to the disc \z\ < 2, we estimate this integral in terms of s. We have 
Kl = \[\Q\\dz\ : and thus the integrand is at most A(s)(\(\ + i)-( d + 4 )/( d + 2 ) 
on any s-admissible curve. Now, an s-admissible curve can be parametrized 
as ( &~ 1 (£ + if{0) '■ ~ oc < £ < °°> where / is a function whose derivative 
satisfies < C(s). So the integral does not exceed 

/CO i 
dei + i)-M/Wy/i + c 2 ( s )de- 

So, the sup of these integrals defining h is bounded by a constant which 
depends only on s, for all Q in a neighborhood of Qd^. 

Proof of Lemma 3. We use the change of the variables which is due to 
Liouville (see, for example, [7]): 

Then the differential equation (THj) becomes 

u>" = W + gu>, (18) 

where 
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where the prime stands for the differentiation with respect to z. 
Now we consider the integral equation 

"(0 = e hC + ^ £je h(C - t] ~ e h ^)g(t)w(t)dt } (19) 

where the path of integration 7 is a part of the <3>-image of an s-admissible 
curve in D. It is easy to verify directly that every analytic solution of this 
integral equation satisfies (TT8]) . (To derive this integral equation one "solves" 
( TT8|) by the method of variation of constants, considering gw as a given 
function, see [7]). 

Putting in (TT9l w = e h ^W, we obtain 

W(C) = 1 + i £jl - e 2h ^)g{t)W{t)dt. (20) 

This we solve by the method of successive approximation: set Wq = and 
define W n +i = 1 + F(W n ), where F is the integral operator in the right hand 
side of ( 1201 . Then we have for the sup- norms on 7: 

\\W n+1 - W n \\, < 4max(l + \e 2h ^\) \\W n - f \g(t)\ \dt\. 

Property (ii) of admissible curves implies that 7 intersects every vertical line 
at most once, so Re (t — () < on 7, and 

| 1 + e 2h(t-f)| < 2 , 

Thus for h > h , the W n 's converge geometrically, uniformly on $(K) to a 
solution W of the integral equation fTSOj) . We have 

W = Wx + (W 2 - W x ) + (W 3 - W 2 ) + . . . = 1 + e, 

where e satisfies (fl6j) . This proves the lemma. 

Completion of the proof of Theorem 2. It remains to prove the statement 
about the uniform convergence. As admissible subsets of canonical regions 
cover Q, it is enough to prove that the uniform convergence holds on each 
admissible subset. Let K be an s-admissible subset of a canonical region D 
containing u + or u~. Our functions Y n satisfy differential equations ([7j). For 
every sufficiently large n, we set h n — \k n \ and 

Q* n (z) = exp(2iargA ; „)(P(A^) - A„), 
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so that 

K = KQlYn- (21) 

In view of (J7J) and aigk n — > 0, we conclude that Q* n — > Q^i as n — > oo 
and all conditions of Lemma 2 are satisfied. When n is large enough, we 
apply Lemma 2 to Q* n and conclude that there exists a canonical region for 
Q* n dz 2 which contains K. According to the Remark after Lemma 3, there is 
a bound for h in Lemma 3 that is independent of n. Applying Lemma 3 
to this canonical region and equation fl2"TT) . we obtain a solution y* of the 
equation (12"TT) which satisfies ffTBT) with h = h n and some branch $ n of 

y \fo*Mdt. 

In particular, this solution ?/* tends to zero on the anti-Stokes direction in 
lo + or u~ . Thus it is proportional to the eigenfunction. Our normalization 
of the eigenfunction implies that the coefficient of proportionality satisfies 
log \ c n \ = 0{n). So we obtain that 

-!-log \Y n \ -> u(z) + c, 
K 

uniformly on K, where c may depend on K. As the limit function exists and 
is continuous in the whole plane, and u is single- valued, we conclude that 
this last relation holds in the whole plane in the sense of D' with the same 
constant c. Comparing the values near we conclude that c = 0. 
This completes the proof. 

Remark. When we were completing this paper, a preprint of Hezari was 
published [8] where similar methods are used to study the self-adjoint eigen- 
value problem 

y" = A 2 g(x)?/, 

with a real polynomial q, q(x) — > +oo as \x\ — > oo and with the usual bound- 
ary conditions on the real axis. Hezari shows that possible limit distributions 
of zeros of eigenf unctions y n are supported by some Stokes lines similar to 
our "exceptional set" . However, in his setting, a single limit distribution as 
n — ^ might not exist, as he shows by an example. 

The main reason why the limit distribution exists in our case is Theo- 
rem 1, which describes the topology of a very special Stokes complex, while 
in Hezari's work an arbitrary Stokes complex symmetric with respect to the 
real line is involved. 
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